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A  conditionally  distribution-free  test  is  proposed  for  testing  the 
symmetry  of  a  bivariate  distribution  function  with  observations  which  are 
subject  to  arbitrary  right  censorship.  In  a  numerical  study,  this  new 
test  is  shown  to  be  more  powerful  than  the  sign  test  under  Marshall  and 
Olkin's  (1967)  bivariate  exponential  model. 


Key  words  and  phrases:  Right  censorship;  Marshall  and  Olkln  bivariate 
exponential  model;  permutation  test. 


1.  INTRODUCTION 


In  many  situations,  the  comparisons  between  two  treatments  based  on 
paired  observations  which  are  censored  In  one  or  both  variates  may  arise. 

For  example,  Hammond  (1964)  has  used  a  matched  pair  analysis  to  study  smoking 
in  relation  to  mortality  in  the  United  States.  Batchelor  and  Hackett  (1970) 
gave  a  comparison  of  the  survival  times  between  HL-A  closely  matched  and 
poorly  matched  skin  allografts  on  the  same  badly  burned  patients.  Also,  in 
life  testing  it  may  be  desirable  to  compare  the  life  times  of  two  components 
in  a  system. 

Suppose  that 


0.2°.y2V 


«°.yV 

n  n 


(1.1) 


are  independent,  identically  distributed  random  vectors,  having  H^(s,t)  as 
their  distribution  function  (d.f.)  and  having  F^(s)  and  G^(t)  as  their 
marginal  d.f.'s,  respectively,  where  '  denotes  vector  transpose.  The  null 
hypothesis,  which  is  to  be  tested  is 


H  :  H°(s,t)  -  H^(t,s)  ,  for  (s,t)  e  R2  . 

o 

Since  and  Y^°  may  be  censored  from  the  right  by  variables  and 

Vj  ,  respectively,  (1.1)  cannot  always  be  observed.  The  observations  avail¬ 
able  to  the  experimenter  actually  consist  of  the  minima 


Xj  -  min(X]w,U1)  ,  •••  ,  Xn 


min(Y.  ,V.)  ,  •*•  ,  Y 
ii  n 


min(X,U) 
n  n 


min(Y\V  )  , 
n  n 


(1.2) 


and 


two  random  sequences  {6,,  ,6  }  and  {c, •  ***  ,£  }  ,  where 

in  in 


2 


0 

1,  if  XA  -  x”  , 
.  o  .  if  X4  <  X®  , 

1,  if  Tj-f®  . 

0 ,  if  yi*. 


For  i  j*  J  ,  the  censoring  variables  and  are  assumed  to  be  inde¬ 
pendent  random  variables  with  a  common  d.f.  J  .  To  have  the  same  censoring 
mechanism  for  both  variates  is  quite  common  in  paired  studies.  It  is  also 

i  0  0' 

assumed  that  (U^.V^  and  (X^  ,Y^  )  are  independent,  1  -  1,2,  •••  ,n  . 

In  the  parametric  case,  the  procedure  for  testing  is  rather  compli¬ 

cated  and  no  useful  results  have  been  obtained.  However,  Holt  and  Prentice 
(1974)  used  the  proportional  hazards  model  (Cox  (1972))  to  analyse  the  data 
by  Batchelor  and  Hackett  (1970),  and  Wei  (1979)  proposed  an  asymptotically 
distribution-free  test  for  Hq  based  on  paired  observations  which  are  subject 
to  arbitrary  right  censorship.  Since  the  sample  size  in  paired  studies  is 
frequently  small,  a  distribution-free  test  is  highly  desirable.  Although 
the  sign  test  is  a  conditionally  distribution-free  test  for  testing  H>  ,  it 
is  rather  inefficient  when  there  are  too  many  censored  pairs  in  the  data. 

As  an  extreme  case,  for  the  data  (3+,4)  ,  (6,5+),  (2+,4),  (9,7+),  (8+ ,6+), 
where  'V  denotes  censoring,  the  sign  test  leads  to  no  conclusion  about  the 
null  hypothesis  Ho  . 

In  this  article,  a  conditionally  distribution-free  test  for  Hq  is 
presented  in  Section  2.  In  a  numerical  study,  it  is  shown  that  the  new  test 
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is  more  powerful  than  the  sign  test  under  Marshall  and  Olkin's  bivariate 
exponential  model  (Marshall  and  Olkin  (1967)). 

We  note  that  all  the  results  of  this  article  can  be  easily  extended  to 
the  case  of  arbitrarily  restricted  observations  (Mantel  (1967)). 


2.  THE  TEST  STATISTIC 

Let  Z±  *  X±  and  Zn+i  *  ,  i  *  1,  •  ••  ,n  .  Vie  will  say  that  Z ^ 

is  definitely  greater  than  Z^  if  Z^  >  Z^  and  Z is  observed,  and  Z^ 

is  definitely  less  than  Z_.  if  Z ^  <  Z^  and  Z i  is  observed.  Now,  let 

£^(n^)  ’  i“l»2»*"»n,  be  the  number  of  the  remaining  (2n-l)  Z's  than 

which  X^(Y^)  is  definitely  greater  minus  the  number  than  which  it  is 

definitely  less.  The  original  observations  (1.2)  are  then  replaced  by 
*  > 

(^l,Th.)  »  ***  *(£n»Hn)  *  Under  and  the  assumption  of  an  equal  censoring 

mechanism  for  both  variates,  all  the  arrangements  of  the  form 


(VR2> 


are  equally  likely,  where  (R2i-i»R2i^  “  °r  (ni»C1>  »  i  -  1,  •••  ,n 

The  statistic  proposed  here  for  testing  H  is 

o 

w  -  I  R2W  • 

i-1  1 


Small  or  large  values  of  7  £ .  lead  to  the  rejection  of  H  .  Note  that 

i-1  1  ° 

U 

7  £ ,  is  Gehan's  (1965)  two-sample  statistic.  Mantel  (1967)  gave  a  simple 

i-1 

routine  to  calculate  and  .  It  should  also  be  noted  that  scores  other 


than  Gehan's  can  be  utilized 
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The  drawback  to  any  permutation  test  such  as  W  Is  the  usually  long  and 
tedious  calculations  required  when  the  sample  size  n  is  large.  Fortunately, 
an  asymptotically  distribution-free  test  is  available  for  large  sample  cases 
(Wei  (1979)).  In  the  rest  of  this  article,  we  concentrate  on  the  small- 
sample  performance  of  the  W  test. 


3.  THE  POWER  STUDY 

In  this  section,  we  study  a  special  alternative  hypothesis  : 

F°(s)  *  G°(s)  for  all  s  and  F°(s')  <  G°(s')  for  some  s'  ,  and  compare 
the  W  test  with  the  sign  test  under  Marshall  and  Olkln's  bivariate  expo¬ 
nential  model.  The  survival  function  of  this  model  with  parameters  X^,  X^ , 
and  X^2  can  be  written  as: 

P{X°is,Y°it}  -  exp[-X1s-X2t-X12  max(s,t)]  .  (3.1) 

The  two  marginal  means  are 

EX°  -  •  l/dj  +  Xjj)  and  EY°  -  1/(X2  +  X12)  . 


Under  this  model,  the  hypotheses  to  be  tested  become  :  X^  ■  X2  against 
H1  !  Ai  <  X2  • 

Three  censoring  schemes  are  considered  in  this  comparison: 

(A)  J(s)  is  a  uniform  distribution  over  (0,EX^)  ; 

(B)  J(s)  is  a  uniform  distribution  over  (0,2EX^)  ;  and 

(C)  J(s)  is  a  uniform  distribution  over  (0,4EX®)  . 

In  this  numerical  study,  the  censoring  variables  and  are  assumed  to 
be  Independent. 
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For  the  sign  test  and  the  W  test,  the  proportions  of  times  in  the 

1,000  Monte  Carlo  samples  generated  that  H  was  rejected  at  the  a  -  .05 

o 

level  were  calculated  for  samples  of  sizes  n  *  10  and  15  from  (3.1)  with 
various  values  of  X^,  X^,  and  X^  •  Tables  1  and  2  give  the  results.  As 
we  expected,  the  W  test  is  uniformly  more  powerful  than  the  sign  test.  In 
addition  to  the  drawback  which  was  illustrated  by  an  example  in  Section  1, 
another  disadvantage  of  the  sign  test  is  the  actual  probability  of  Type  I 
error  is  far  below  the  specified  a  value.  For  example,  when  n  ■  10  , 
under  the  severe  censorship  (A),  the  empirical  levels  of  the  sign  test  are 
only  .01  as  compared  with  the  nominal  value  a  «  .05  . 
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